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Energy variation of the density of states (DOS) has been calculated in the superconduc- 
tor/ferromagnet/ferromagnet /superconductor structure (SFFS) in the frame of Gorkov equations 
taking into account the s-d electron scattering in the ferromagnetic layers. DOS behavior is pre- 
sented for the antiparallel and parallel magnetic moments alignment of two adjacent F layers. The 
cases of small and large values of exchange ferromagnetic field are discussed. 
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I. INTRODUCTION 

The possibility of the potential application of multi- 
layers with the superconducting junctions stimulate an 
interest to the study of such a structures. The effect of 
switching between "0" and "7r" states in superconduc- 
tor/ferromagnet/ferromagnet /superconductor (SFFS) 
structures makes them useful for the microelectronic 
engineering. Practically it is possible to change the 
mutual orientation of the thin ferromagnetic layers 
magnetization by applying the external magnetic field. 
Two ferromagnetic layers separated by a nonmagnetic 
spacer FNF are known as " spin- valve" . The combined 
SFFS structure with superconducting electrodes gives 
more interesting results. In the parallel configuration for 
the some ferromagnetic thickness values the Josephson 
current takes a zero value. At the same time for the 
antiparallel alignment of magnetizations the Josephson 
current has a non zero value. Thus, by selecting a 
necessary ferromagnetic layer thickness we can govern 
the superconducting current flow. 

The presence of two ferromagnetic layers with the di- 
rigible magnetization alignment let study two opposite 
cases which have their peculiarities. For the parallel mag- 
netization alignment it is the reversion of the Josephson 
current sign (n state junction) that has been theoretically 
described Q and experimentally observed 0, Q , and for 
the antiparallel case the possible current divergence. 

Several authors involved in the studying of SFIFS 
structures (I is the insulating barier) touched upon 
a question of the possible enhancement of the criti- 
cal Josephson current for the antiparallel magnetization 
alignment of the neighboring F layers under some con- 
ditions. In it was found that in the case of the an- 
tiparallel orientation the critical current increases at low 
temperatures with increasing exchange field h and at zero 
temperature has a singularity when h equals to the su- 
perconducting energy gap A. Authors of also pre- 
dict the critical current enhancement for the antiparal- 
lel alignment in some interval of exchange field values. 
The possibility of the critical current enhancement by 
the exchange field in SFIFS junctions with thin F layers 
with antiparallel magnetization directions was discussed 



in the regimes os small S layer thickness |j] and bulk 
S electrodes 0. Golubov et al. Q assumed that the 
case of thin S layers 3 was studied in the frame of an 
idealized model in the tunnelling limit, which leads to a 
divergency of the critical current at the zero temperature, 
and for the bulk S case || an approximate method was 
used, so that a part of the results was obtained beyond 
its applicability range. Golubov et al. explain the cur- 
rent enhancement by the singularity in density of states 
(DOS) which is shifted to the Fermi energy level. It hap- 
pens when the effective energy shift in the ferromagnets 
due to the exchange field becomes equal to a local value 
of the effective energy gap induced into F layers. They 
note that in the models with DOS of the BCS type this 
leads to a logarithmic divergency of the critical current 
in the antiparallel case at zero temperature. 

First of all, it should be emphasized that in most the- 
oretical papers the SFFS structures were studied in the 
so-called dirty limit and for low energies, that is using 
the quasiclassical Usadel equations in the context of one- 
band model of ferromagnetic metal. Such approach docs 
not take into account the s-d scattering of electrons which 
is the main scattering mechanism responsible for the ki- 
netic properties of 3d-metals and their conductivity. Ev- 
idently, the limitation connected with the Usadel equa- 
tions may be avoided by solving the full Gorkov equations 
which in addition take into account the s-d scattering of 
conducting electrons. 

In the paper Q the Josephson current in SFFS struc- 
ture had been calculated in the frame of Gorkov equa- 
tions taking into account the two-band model of a ferro- 
magnet with conducting s electrons and almost localized 
d electrons. S layers were assumed as a simple s-wave 
superconductors. For the antiparallel configuration the 
Josephson current has no singularity at zero tempera- 
ture. Consequently, from the physical point of view the 
results derived in jj, 0, 0] stay unclear and this question 
requires a more careful study. For that reason in this 
paper we calculate and discuss DOS for both (parallel 
and antiparallel) cases in the frame of Gorkov equations 
taking into account the s-d scattering of electrons. As it 
will be shown DOS does not demonstrate any divergence 
that makes us conclude that s-d scattering effectively de- 
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stroys the BCS correlation, hence, the Josephson current 
enhancement is suppressed. 

II. GORKOV EQUATIONS 

We consider the SFFS structure which consists of two 
semi-infinite superconducting electrodes S separated by 
two thin ferromagnetic layers F with the interfaces per- 
pendicular to the z axis. Although the Green function 
method can be applied in more general case, for simplic- 



ity we assume the F and S materials to be equivalent 
on both sides of structure, and both F layers have the 
same thickness a which is supposed to be much smaller 
than the in-plane dimension of the structure. The S/F 
interfaces are placed at z = ±a. In the mixed ( k , z) rep- 
resentation ( k is the quasi-momentum component in the 
XY plane of layers) the Gorkov equations for the normal 
G|| and anomalous F" Green functions have the follow- 
ing form: 
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where lo = 7rT(2n+l) is the Matsubara frequency at tem- 
perature T, e F {z) = e s {6(-a - z) + 9(z - a)} + e f 9(z + 
a)0(a — z), e s ^' the Fermi energy of the superconduc- 
tor (ferromagnet) , 6(z) the step function (9(z) = 1, if 
z > 0, and 8(z) = 0, if z < 0), h = k B = 1. We con- 
sider that the Cooper pairing is absent in the F layer 
and the superconducting order parameter A is equal to 
zero in the ferromagnet, A(z) = Aexp(itpx)8(— a — z) + 
Aexp(iif2)0(z + a) (^1,2 is the order parameter phase 
in the left and right superconducting electrodes, corre- 
spondingly). The exchange magnetic field in the F layer: 
h(z) = h9{z + a)6{-z) ± h6(z)0(a -z),h> 0, sign (+) 
corresponds to the parallel alignment of F layer magne- 
tizations and (— ) sign should be taken for the antiparal- 
lel configuration. The s-d scattering is described in the 
framework of the fist Born approximation with impurity 
concentration xq(z) = xq8(z + a) 9 (a — z) and impurity 
potential value 7. G^^^(z, z, lo) are the diagonal Green 
functions of localized d electrons. Here we neglect the s-s 
electron scattering as it is small in comparison with the 
s-d scattering which are the main mechanism destroying 
induced superconductivity in the F layers. 

In this case the Gorkov equations become linear and 
the solutions represent a set of plane waves with wave 
vectors: 



k m = J2m{ef±h-£±iu± ^j, 
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to the the S layer electrons. For the parallel configura- 
tion ki = fc4, ki = h§. Here t-j-q) is the spin up (f) 
and spin down (J,) electron lifetime in the ferromagnet, 
Vnu) = - 2 *o~/ 2 ImG$H l y 

The procedure of DOS calculation reduces to the di- 
agonalization of the imaginary part of the normal Green 
function 
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K(IkG^,^(z, z, iu> — ► uj). (3) 



Because of unhandiness of the expressions for the 
Green functions here we cite as an example the normal 
Green function only of the left F layer for spin up elec- 
trons: 
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where f± and fa are the functions of electron momenta, 
lo and A which are given in the appendix. 

In order to discuss the connection between the behav- 
ior of DOS and Josephson current which served as a mo- 
tivation force of the present study we rewrite here the 
expression for the current derived earlier (see 0) m th° 
frame of Gorkov equation for the antiparallel and paral- 
lel orientation of the F-layer magnetic moments, corre- 
spondingly: 
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where ku2) corresponds to the majority (minority) spin 
electron momentum in the left F layer for the antiparal- 
lel configuration (at the same time in the right F layer 
with the opposite magnetization direction majority and 
minority momenta we label as k^ and k§), k$ corresponds 



MO)^^/^ ./ 1 -^. , (6) 

where if is the phase difference between right and left 
superconducting electrodes order parameters, R and r 
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arc reflection coefficients, Jap and fp arc certain oscil- 
latory functions which depend on the ferromagnetic lay- 
ers thickness a, energy parameter lu and superconducting 
energy gap A. There is no singularity in the Josephson 
current for the antiparallel case at low temperatures (see 
Fig. nj predicted previously in 0, |(| . 
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Figure 1: Josephson current J versus F layers thickness a. 
Fermi momenta in F layer fc T = 1 A" 1 , k L = 0.8 A -1 (fc T( - L) = 
\j2m(E$ ± h); in S layer k s = 1.2 A -1 . Free path of electrons 
in F layer if = 100 A and l± — 60 A. The critical temperature 
of the superconducting metal is T c — 4 K. 



III. ANTIPARALLEL CASE 

Antiparallel case. In Ref.[(J the authors derive the ex- 
pression for DOS which demonstrate the energy renor- 
malization due to exchange field in SFIFS (I - insulating 
barier). For a certain exchange field value h DOS expres- 
sion yields the singularity which is shifted to the Fermi 
level that leads to the logarithmic divergency of the crit- 
ical current at low temperatures. However, the authors 
avoid the singularity problem by taking the low barrier 
transparency which makes peaks broader and in such a 
way suppresses the singularity. 

In Ref.|| energy variation of DOS was investigated 
in the S/F structure. This case resembles more to the 
SFFS structure with parallel magnetizations of F layers, 
but here we can note that just as in Ref.[(| DOS has a 
peak and a region of zero value within the energy gap 
in much the same way as in S/N structures ( N is the 
normal metal) 0- 

Both calculations were made for the small value of the 
exchange field h ~ A. Our computation for the same 
range of exchange field is depicted in Fig. [3 As one can 
see the DOS has a BCS behavior outside the energy gap 
and usual divergence at uj — A. However there in no any 
singularities or peaks within the energy gap and DOS 
never equals to zero. Last peculiarity may be explained 



by the destructive role of s-d electron scattering in F lay- 
ers which destroys the Cooper pairs. The curves a,b,c 
were plotted for the different values of the free path i-m) 
of spin up(down) electrons in the F layer, l^^ = vft^^, 
vf is the Fermi velocity in the F layer. For the small 
values of free path (i.e. in the case of the strong s-d scat- 
tering) the hollow becomes smaller and DOS approaches 
the constant as in the bulk ferromagnet. The inset of Fig. 
|21depicts energy DOS variation within the energy gap for 
the different F layers thickness a. For the decreasing val- 
ues of a DOS becomes more close to zero as in a bulk 
superconductor. The peaks in DOS appear for the large 




Figure 2: Energy variation of DOS at the left F/S border (z = 
—a) for the antiparallel case for the small value of exchange 
energy h ~ A, A = 1.4 10~ 3 eF. a) a = 7 A, Z T = 500 A; b) 
a = 7 A, ly = 100 A; c) a = 7 A, ly = 10 A; d) a = 15 A, 
l } = 500 A; e) a = 25 A, Z T = 500 A. 
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Figure 3: Energy variation of DOS at the F/S border (z — 
a) for the antiparallel alignment of F-layers magnetization. 
Fermi momentum of electrons in S-layer: k s = 1 A -1 . In the 
F-layer for the electrons with spin up: k} = 1 A -1 , k^^ = 
y/2m(el ± h). a) Z T = 500 A; b) Z T = 100 A. 

values of the ferromagnetic exchange field h ~ 10 3 A (the 
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case of a strong fcrromagnct). Outside the energy gap 
DOS stays of a BCS type and within the gap DOS has 
two symmetrical peaks (Fig. As in Ref. they 

can be shifted closer to the Fermi energy {u> = 0) but at 
the same time their size diminishes and there is no any 
singularity at the very zero. For the smaller values of free 
path and large values of s-d scattering parameter (Fig. 
Ob) the peaks lose their pointed shape. 

In both cases of small and large values of ferromagnetic 
exchange field DOS does not demonstrate any divergence 
at the Fermi level which is in a good agreement with the 
calculated earlier Josephson current dependencies for the 
antiparallcl alignment of F layers magnetic moments. 

IV. PARALLEL CASE 

For the small values of the exchange field h ~ A in the 
parallel case there are no any significant changes in DOS 
in comparison with the antiparallel case presented in Fig. 
12 But for the large values of h in contrast to the an- 
tiparallcl case multiple electron reflection takes place that 
graphically manifests itself in multiple resonance peaks in 
DOS within the energy gap (see Fig. 0}. 
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Figure 4: Energy variation of DOS at z — for the paral- 
lel alignment of F-layers magnetization. Fermi momentum of 
electrons in S-layer: k s = 1 A -1 . In the F- layer for the elec- 
trons with spin up: fc T = 1 A -1 and exchange field h ~ 10 3 A. 
a) "0" state: a = 14 A, Z T = 500 A; b) "0" state: a = 14 A, 
Z T = 100 A; c) "tt" state: a = 19 A, Z T = 500 A; d) "tt" state: 
a = 19 A, Z T = 100 A. Dashed line: DOS for spin up Nf. 
Dotted line: DOS for spin down iV;. Solid line: JV T + N x . 

This fact may be explained in the following way. The 
antiparallcl alignment is more favorable for the relatively 
free diffusion of the Cooper pair electrons which have the 
opposite spin direction. Possible destruction of the pair 



in the first F layer compensates by the presence of the 
second F layer with the opposite magnetic moment align- 
ment and reversed majority and minority spins. While in 
the parallel case for the large values of the exchange field 
/i > A, when spin up and spin down electron energies 
differ notably only one electron in the pair whose Fermi 
energy is the same in the F and S layers may travel rel- 
atively free. Hence, the initial pair destroys rapidly and 
the multiple reflections becomes possible as it happens in 
the potential well. 

DOS dependencies in Fig. 01 are plotted for the differ- 
ent values of free path. For increasing value of s-d scat- 
tering in the F layers multiple irregularities arc smoothed 
over and energy DOS variations resemble to those calcu- 
lated near the S /F border of bilayer in m the frame 
of Gorkov equations. So called "0" and "7r" states are 
presented in Ref. 0. Due to the oscillations of the su- 
perconducting order parameter induced in the F layers in 
the "7r" state the DOS shape is reversed with the respect 
to the normal "0" state. Instead of usual BCS energy 
dependence with two peaks at u> = ±A DOS has two 
dips at these energies values. As at was discussed earlier 
the Josephson current varies its sign in the parallel case 
[ME10]|< (+) corresponds to the "0" phase and (— ) to the 
"tt" phase. We have compared DOS behavior with pre- 
sented in Fig. ^ Josephson current oscillation from " 0" 
to "7r" state in the parallel case. Thick points in Fig. 
correspond to the F layers thicknesses a which had been 
chosen for that comparison. Indeed, for these values of a 
DOS demonstrates the same transition between "0" and 
"tt" states (Fig. 0]a,b and Fig. E|c,d, correspondingly). 



V. SUMMARY 

In this paper the energy DOS variation had been stud- 
ied in the SFFS structure for the antiparallel and parallel 
magnetic moments alignment of two adjacent F layers. 
It was shown that there in no any singularity in DOS 
in the antiparallel case that is in a good agreement with 
calculated earlier Josephson current behavior in the same 
model. In the parallel case the presence of "0" and "tt" 
states had been shown. 
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APPENDIX: EXPRESSIONS FOR THE 
COEFFICIENTS IN EQ. © 
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fi,2 = TT-H fc i ± fc 4) (Ai(k 3 + ki) + A 2 (k 4 - k 3 )) + 7n-( fc i T k 4 ) {Ai(k 4 - k 3 ) + A 2 (k 4 + k%)) , 
Ak\k 4 \k\k 4 

Ai, 2 = ^ [±(B a *i - B 1 kl)(k 3 t 2 ,i + k 2 t 4 , 3 ) ± A 2 (B 1 fc 3 + S 2 fc 1 )(fc 2 t 4 ,3 - fc|i 2 ,i)± 

±AiA 2 *i(fc 3 + K){B 2 c 2 ,i - 5x04,3)] , 

den = fc 2 (fc3 + fcl) ^2*3 - *i*4) +AiA 2 fci(fc3 + fcg)(ciC4 + c 2 c 3 )+ 

+Ai ((fcic 4 - fc 3 c 2 )(Mi + fc 2 i 3 ) + (fcgCi - kic 3 )(k 3 t 2 + k 2 t 4 )) + 
+A 2 ((fcic 4 + k 3 c 2 )(k 2 t 3 - k 3 h) + (fc 3 ci - fcic 3 )(fc|t 2 - fc 2 i 4 )) , 



Si, 2 = i( e - iz ' kl e- iakl Te iz ' kl e iak ^)/2k u A 1:2 = ie - **' 1 

c i,3 = ^7 [e- wfe4 (fc 3 + fc 4 )(fci+A; 4 ) + e^(fc 4 -fc 3 )(fc 1 -A :4 )]± 
± & [e-^ih + h^h-k^ + e^ih-k^ih+h)] , 

c 2,4 = s ^[e- taki (k 4 -k* 3 ){k 1 +k 4 ) + e tak *(k 4 + k* 3 )(k 1 -k 4 )]± 
±& [e~ wfe4 (fc 4 - fc 3 *)(fci - fc 4 ) + e» fe *(fc 4 + k* 3 )(h + k 4 )} , 

ti, 3 = A 1 |^[e 8a ' £2 (A :2 -A :3 )(fc 2 + fc5) + e- iafc2 (fc2 + fc3)(fc5-fc 2 )]± 
±A X ^^ [e 4afc2 (fc 2 -fc 3 )(fc 5 -fe) + e-^(fc 2 + A ;3 )(fc 2 + fc 5 )] , 

t 2 , 4 = -A 2 |^[e Mfe2 (fc 2 + fc3*)(fc2 + fc 5 ) + e- ia ^(fc2-^)(fc5-A;2)]T 
T A 2 ^^ [e"^(fc 2 + fc 3 *)(fc 5 - fc 2 ) + e-- fe *(fc 2 - k* 3 )(k 2 + k 5 )]. 
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